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Appendix III 

III.l Strain Energy Release Rate 

In this analysis, a delamination between belt and core sections is assumed to grow 
parallel to the belt direction in the tapered and uniform sections. These delaimna- 
tions in each section are denoted by a and 6 respectively. The core section in the 
taper portion is modelled by two equivalent sublaminates. The stiffness properties 
are smeared to obtain the effective cracked and uncracked stiffnesses which are des- 
ignated by Ay and Ac as shown in Figure III.l. These stiffnesses change from one 
plj' drop group to another with crack growth o by experiencing a sudden change 
at discrete locations. Therefore Ay and A<- can be represented in three consecutive 
regions as follows, 


• Region 1; 0 < a < / 


Ay 


d ol — a 

d I I I I I q 
^BD ^2 ^3 


(HI.l) 


Ac = As 


• Region 2: I < a < 21 


Ay 


d + 31 — a 

d I / 4 2/-g 

Abd Ai Aj 


. a -f- 6 

A2 ^ A^ 
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• Region 3: 21 < a < Zl 


Ac 


d Zl — o, 

ZZTTX 

Abd Ai 


a + b 


a-1l 

At 




l±k 


(IIL5) 

(111.6) 


where 

h = ply thickness 

d = length of uniform thick portion . 

/ = distance between two consecutive ply drop locations 

.4i = GhQ^'- ^ 2hQ° 

A 2 = 4hQ^'^ + 2/iQ® 

A3 = + 2hQ° 

Abd- 7kQ^ + 2hQ^^ 

Q° = Qii of a 0 degree ply 
Q'*® = Q\\ of a ±45 degree ply 

Geometry of the sublaminate model is shown in Figure (III-l) 
Also axial stiffnesses Ajj, A, and Ap are given by 

d ± Zl — <z 

Ab = 


d , 3/-0 
Abd ^bt 


(III.7) 
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Figure III.l: Geometry of the Sublaminate Model 
Ap = .43 

- 4 , = Abd 

w^here 

Abt =Taper belt stiffness 

For a membrane behavior, equilibrium equations are reduced to 

i\T,,= 0 

and the displacement field is assumed to be 

u(i,2) = C7(i) 


(III.S) 

(III.9) 


(111.10) 

(111.11) 


w — 0 


(III.12) 
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The constitutive relations are represented by 


N = AnU^ (IIL13) 

The stress and displacement fields, are determined based on the stiffnesses derived 
in Equations(IILl-III.9). In this model, load is shared by the core and the belt 
portions according to their respective stiffness ratios 


Pi = 


PAb 
Ab + Au 


(III.14) 


P2 


PA, 

Ab + Au 


(III.15) 


where P is half of the total axial load applied at the ends. 

Using the Equations (III. 10), (III. 13), and the expressions for Pj and P 2 from 
Equations (III. 14). (III. 15) the axial displacement at x = c can be written as 


PAbc P{d -T 31 -r b) / Ab Ab \ 

A.(Ab + Au) ^ (Ab + Au) VA^ ~ 17/ 

P{d + 3l-a) / _ Ab \ 

{Ab A u) \ Ab\) 


(III.16) 


P A,c P{d + 3/ + 6) / Au Au \ 
Aa(AB + Au) (Ab + Au) VA,. Af) 

Pid + 3l- a) / _ 

(Ab+Au) I aJ 


(III.17) 


where Abi is the belt stiffness in the pop-off region as shown in Figure Ill.l. 

A three-dimensional transformation is required in order to estimate the effective 
ZLxial stiffness of the belt region Ab and .4 bi- This is due to the belt layup and 
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the orientation of the different belt portions to the loading axis as shown in Figure 
III.l. The three-dimensional transformation is presented in section III.3. 

The tapered lauminate is assumed to be fixed at i = 0. Therefore the external 
work done is given by 


W = + P2Ue 


(III.18) 


Substitute from Equations (III.14) through (III. 17) into Equation (III. 18) to get 


W 


(Ab + A^Y 


{c — d — Zl — b) ~A^^ {d A Zl — cl){Ab + .-4V,) 

(III.19) 


/4i '' 


The strain energ}' release rate G due to the external work done is determined 


bv 


G = 


1 dW 
2P2 dA 


(III.20) 


where A is the delamination surface area. G is calculated for delamination lengths 
ranging from 0 to 60h. In the analysis, S2/SP250 Glass- Epoxy is used. Its properties 
are given in Table III.l. 


Table III.l: Material Properties of S2/SP250 Glass-Epoxy 


Ell (MSI) 

E 22 (MSI) 

Gi2 (MSI) 

Gi3 (MSI) 

G 23 (MSI) 

1^12 

7.3 

2.1 

0.87 

0.5 

1 

] 

0.5 

0.275 
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II1.2 Interlaminar Stresses 

In this part, an analysis for the interlaminar stresses in the belt-core interface in 
the tapered section will is developed. 

The simple analytical model assumes a beam model for the belt in the tapered 
section which is shown in Figure III.2 . Material and geometric discontinuities are 
modelled as extensional fc,; and concentrated shear springs (i=l-4) as shown in 
Figure III.3. The resin pockets are assumed to be subjected primarily to shear stress 
and they are represented by a distributed shear spring with a constant stiffness G. 
The effect of the core is incorporated as elastic supports on the beam-belt model. 

A minimimi complementary potential energy formulation is used to estimate 
the interlaminar stresses. The total complementary potential energy consists of 
bending, shear and extensional energy contributions, 


= Hfc + n, -4- n, + n,. 


(III.21) 


where I!;,, II 4 . !!«, 11^ represent bending, shear and extensional energy components 
and energy stored in elastic springs, respectively. These are given as. 


Hi. 


1 M-{s) 

2 Jo -Dll 


(III.22) 






1 /-3' N-[s) 


ds 


>■11 


(III.23) 


(III.24) 




2 Jo G, 2fc, 


+ + ^ ■ IL + IL + H.IL 

2^2 2 1^3 2k ^ 2gi 2 q 2 2g^ 2g^ 


(III.25) 




Figure III.2: Geometry of the Model 
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where Ri, T{ (i=l,2,3,4) are unknowns. The constant shear stress, c, due to resin 
filler is an additional unknown. The total number of unknowns in this formulation 
is nine. These unknowns axe constrained by following equilibrium equations. 

Ri = -Rz - 2 R^ + 2Nn + N ^2 - - N^i ) (III.26) 

R 2 = 2Rz + ZR 4 - SATij + ^{Nn - ) (III.27) 

Ti = -Ti -Tz-Ti- 3c/ - ^^2 -f (III.2S) 

where A/n, Nn, N 21 and N 22 denote the components of the extensional load at two 
ends of the belt section. 

The bending moment, shear force and axial force in each of the three ply drop 
regions are written as 

• Region 1: 0 < s < / 

M(^s) = —N\ 2 ^ + + R 4 S + T4~ (III. 29) 

V(s) = Ni 2 - Ra (III.30) 


N{s) = Nu -cs-Ta 


(III.31) 
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• Region 2: I < s < 21 


M(s) = -N^ 2 S + JS + {Rz + Ra)s - Rzl + {Tz + T^)^- 


V{s) = Nr 2 -Rz-R4 


N{s) = Nn -cs-n-T^ 


• Region 3: 21 < s < Zl 


M(s) = -Ni,s + -f (-R, -hRz^ RaU + (2R2 - Rz)/ -f (T 2 + Ts + T,]^ 


V(S) = A*i2 - R2 -R3-R4 


N{S) = JV„ -C 3 -T 2 -T 3 -T^ 


Therefore the bending energy in Equation (III.22) can be written as 


Hfc = 


— f‘\- 

2Z?ii Jo 




t2 


ds 


+ ^J, [--^12" + f ^ + (^3 + Ra)s- Rzl + {Tz + T,)^- 

1 f .. ci 

^w7,L + + 


-R, - 2R, + N„ - -(A’„ - A', 


(111.32) 

(111.33) 

(111.34) 

(111.35) 

(111.36) 

(111.37) 


s 


11 


3i?3 + QRa — QNi2 + ^(-^11 ~ -^21 )j I 

+(^2 + ?3 + r^) - 1 ds 


Similarly for the shear energy 


(III.38) 


n, — f" {Ni 2 — R4)^ ds + r {Ni 2 — Rs — a 


ds (III.39) 


where 


_ _3 


SG-,A • 


The energy of extensioiial loads can be expressed by 


n, = 


1 r' 


2 A 


11 


jT ds + ^ j'\Nn -CS-T 2 - T 4 )- 


ds 


+ -CS-T2-T,- T^f ds 

2 ,Aii J 21 

The energy stored in the elastic springs is written as 


TT 3 1 

2G’2^2fciL 


R3 — 2R4 + 2Ni2 + N22 — ^(-^11 ~ -^21 ) 


+ 21^^ + 3R4 — ^ Mi + j-^^2 ~ 3 A’^i2 + ~ •^'^ 21 ) 


(III.40) 
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Ti Ti Tl 

2^2 2^3 2^4 


(III.41) 


The complementary potential energy in Equations (III.38) through (III.41) is ex- 
pressed in terms of 6 unknowns, namely ijg, Ti (i=2,3,4) and c. By minimizing 
these expressions the following system of linear equations is obtained 


/27t2/3 ^ 9^\ 

1,12-D 


„ 5<P „ 


SD 





Ts 


'9tH^ 
, SD 


5tP 

2D 


+ -\T, = ^ W. -N,,)+ ■■^(Nn - N,,) 


tH\ 

3d' 


(III.42) 


iP (2P 1 4 1 \ 

o ^ I TT: + T~'j-7“ + — l-R3-r 


D 


ZD k-i A’2 k-z 




.2£> ^•l kz^ 




tl^ 


tP 


ZP 


j rp , ZI—T — I ^ ^ 1 A- 2 i \^22 


+ 


f iP t t \ , 
\12T> 2k^l ~ kzl j 


(III.43) 


btP f ZP 2 6\ /4P 4 9 1 \ tp ^ii- 

' 2 B (25 + 1 ;) ("d + + 1 ; + iti) ^ * Id 


3t/2 

21 ) 


_ f4P 4 9 \ 27 V 22 


(III.44) 
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'btH^ Zl\ tP 


( tn 1 \ „ tp 1 iH 

■" iis + j ^•' = 55*" + - *=>) + 5 d(*" - *■«) (Hi-«) 


'PP Zl\ tP ^ 5tp (PI l\ (PI 1 1 


AD 


(Pi l \ ^ 5tP , 1 Pi 

\ 2 D gj ~ Id - -^21 ) + ^(^u - ^'21 ) (III. 46 ) 


9PP Zi' 




3f./2 


Pi 1 


/ /-/ 1 


SD ^ly ' ■ 2D"'’ 2Z>-- V4I> 9x1'^ \2D gj 
( Zt^l 1 1 \ _ Zip 1 

Pi 


+ 


16 D 


(^n-iV2i) 


(III.47) 


Tlie concentrated, norinal and shear forces at the ply drop regions and the inter- 
laminar shear in the resin filler are estimated b}^ solving the simultaneous system of 
equations in (III.42) through (III.47) and using Equation (IIL26) through (III.2S). 
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111.3 3-D Transformation of Stiffnesses 

It has been determined that a three dimensional transformation of stiffnesses is 
required in order to estimate the effective axial stiffness of the belt regions, Ab aiid 
Ab\- This is due to the belt layup and the orientation of the different belt portions 
to the loading axis as shown in Figure III.4. 

The loading axis corresponds to axis 1 in the 123 coordinate system which is the 
transformed system. The principal material coordinates are denoted by l',2’ and 
3'. 

The stress-strain relationships in the principal materiad coordinates for an or- 


thotropic laminate 

axe given by 



= (QJsxe {€}exi 

(III.48) 

where 



On 

= (1 — ?-'23^'3l)l jE'13 

(111.49) 

Q 22 

= (1 - E 22 

(III.50) 

Qz3 

— (1 — 1 ^ 121 ^ 21 ) 1 ' E33 

(III..51) 

Qi2 

= ( 1^21 + ^23^^31 )F-^n = (^12 + ^ 13 ^ 32 )^ ■£'22 

(III.52) 

Qi3 

= ( 1^31 + ^ 21 ^ 32 )V Ell — ( 1^13 + ^ 23 ^ 12 )V E33 

(111-53) 

Q 23 

— (l''32 + 1^12*^31 )1'^ ■£'22 = (*-'23 + 1^21*^13)1^ -£'33 

(III-54) 

Q 44 

= O 23 

(III-55) 

Q 55 

= O 31 

(111-56) 

Qee 

= G'i2 

(III-57) 
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Figure III.4: 

^ = (1 “ — '^23'-^32 — J''31»'l3 — 2 j^12/^23'^31 )'^ (III. 5$) 

The presence of airgle plies in the belt region making an angle B in the l '2 -plaiie 


results in the following constitutive relationship 


K } = m {^'} 

(III.59) 

where the transformed reduced stiffnesses Q^j are given in terms of reduced stiff- 
nesses Qij as 

Qij = c^^Qu + 2c^s^Qi2 + s^Q22 + 4c^5^Qe6 

(III.60) 

Q 22 “ + 2 (^s^Qi 2 + <^Q 22 + 

(III.61) 

Q^2 = + s^)Q^2 + c^^^Qi2 — ^c^s^Qee • 

(III.62) 


^66 — 4c"S^Qll — 8c^S^Qi2 + 4c^S"Q22 + 4(c^ — 5^)^Q66 


(III.63) 
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Qsz 

= Q33 

(III.64) 

Qiz 

= C^QlZ + ^^^23 

(III.65) 

Qzz 

= <S^Qi 3 -t- C^(?23 

(III.66) 

c 

= COS0 


$ 

= sinO 



Any ply in the belt portion of the taper makes an angle /? with the loading axis 
if it is in the uncracked belt portion and an angle a if it is in the cracked belt 
portion. By performing a rotation about the 2-axis, the stiffness along the loading 
axis , takes the form 


1 




(111.67) 


where cr,j and are in 123-axis system and G,, represent the elements of trans- 
formed stiffness matrix in this coordinate system. 

Since we have assumed 


u{x,=) = U{x) 


(111.68) 


and 


U’ = 0 


(111.69) 


For plane stress condition in 1-3 plane (i.e. <Ti 2 = 0 ; i =1,2,3) stress strain 
relations reduce to 


<^11 — ( Gji — C12/C22) ^11 


(111.70) 
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where 

Cn = + 2c -f- + cVQss (III.71) 

Cu = c^Qi2+s^Q23 (III.72) 

C 22 = Q 22 (III. 73) 

where c and s are cosine and sine of the angle which the cracked and uncracked belt 
portions makes with the loading axis. 

The coefficient of cn in Equation (III.70) represents the transformed axial stiff- 
ness. This value is used in the derivation of As and 





